Crosslinked networks of stiff filaments exliibit negative normal stress 
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Motivated by recent experiments siiowing that stiff biopolymer gels exhibit highly unusual negative 
normal elastic stresses, we develop a computational model for stiff polymer networks subject to large 
strains. In all cases, we find that such networks develop normal stresses that are both negative and 
of magnitude comparable to the corresponding shear stress. We find that these normal stresses 
coincide with other nonlinearities in our networks, and specifically with compressive bucking of the 
individual filaments. Our results suggest that negative normal stresses are a characteristic feature 
of stiff (bio)polymer gels that have been shown to exhibit strong nonlinear elastic properties. 
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Networks of semi-flexible polymers such as those that 
make up the cytoskeleton of plant and animal cells 
have been shown to have rich mechanical and rheolog- 
ical properties. One of the hallmarks of their mechan- 
ics is a highly nonlinear elastic response to stress and 
strain [Ij El [3l IH [3 |6] , including dramatic stiffening un- 
der shear. A striking example of such nonlinearity has 
been the recent demonstration in such systems of highly 
unusual negative normal stresses — e.g., in which a sam- 
ple will tend to contract along an axis perpendicular to 
the direction of shear[7]. Normal stresses, in general, 
are a nonlinear phenomenon, since their sign cannot de- 
pend on the direction of shear, for symmetry reasons. 
But, most materials tend to expand when sheared, as 
has been known at least since the classic experiments 
of Poynting nearly 100 years ago[5], in which he showed 
that elastic rods extend axially when twisted. Another 
familiar example of this is the tendency of granular ma- 
terials to dilate when sheared, as can be seen by the fact 
that wet sand tends to dry out around our feet when 
we walk on the beach. Few materials have been found 
to develop negative normal stresses. Liquid crystalline 
polymers [9|, nanotubes|10] and emulsions [llj are exam- 
ples of such systems, which show rather weakly negative 
normal stresses in a range of applied shear rates. By 
contrast, semi-flexible polymer networks exhibit negative 
normal stresses of comparable magnitude to the observed 
shear stress, and this is observed in the elastic response 
of such networks. 

Here, we simulate networks of stiff rods and show that 
these networks generally exhibit negative normal stresses 
(NNS) comparable to or larger than the shear stresses, 
depending on network structure and the strain applied. 
It has previously been shown theoretically that entropic 
effects can result in such NNS effects [7] . We also find such 
anomalous NNS effects in a purely mechanical (ather- 
mal) model of stiff rods. We study the dependence of 
both shear and normal stresses on both the density and 
individual mechanical properties of the constituent fila- 
ments. Although both entropic and enthalpic networks 
can exhibit NNS, their predicted dependence on network 



FIG. 1: A portion of a typical network strained at 7 ~ 0.16. 
Red rods (color online) are stretched segments of filaments 
while blue ones are compressed. Filaments oriented around 
45° are stretched while the ones at about 135° are compressed 
and buckled. Horizontal boundaries are rigidly displaced ac- 
cording to a macroscopic strain 7 while the sides are con- 
nected by periodical boundary conditions. 



properties such as density are opposite, providing a pos- 
sible experimental way to distinguish the important roles 
of entropic and enthalpic effects. 

We construct model networks as follows. A number 
of straight filaments of fixed length L and random ori- 
entation are deposited on. a,W xW square. Every time 
two filaments cross each other they are linked together 
by a free hinge. Filaments that cross the upper or lower 
boundary are rigidly attached to the boundary and dan- 
gling ends are removed, while periodic boundary condi- 
tions are applied to the left and right edges of the square. 
The deposition stops as soon as the desired concentration 
is reached. We characterize this concentration by the ra- 
tio L/ ic, where ic is the average distance between cross- 
links along a filament. The network is then refined by 
adding midpoints between any pair of consecutive cross- 
links along the same filament. An energy functional is 
associated with the position of the modeled points, both 
cross-links and midpoints, in such a way as to implement 
a discrete worm like chain approximation for every fila- 
ment. For any consecutive pair of points we consider the 
stretching energy ^TYstretch as a function of the distance 



£ and the initial distance Iq via a stretching stiffness /z 
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The bending energy is evaluated for every pair of consec 
utive segments along a filament: 

2 



sn^ 



(2) 



where 69 is the angular deflection of one segment relative 
to the other, £' is the average length of the two segments 
and K is the bending stiffness. The total energy is, then, 
the sum of these energies over the entire network. A 
shear of this network is performed by displacing the up- 
per boundary with respect to the lower one, followed by a 
minimization of the elastic energy over the unconstrained 
internal coordinates (cross-links and midpoints) of the 
network using a nonlinear conjugate gradient technique 
[T^ . Finite deformations are accomplished in a stepwise 
manner. The resulting shear (normal) stresses are then 
determined from the forces parallel (perpendicular) to 
the displaced boundaries. 

We simulated networks of size W = 7 oi filaments with 
length L — I, 2 and different network densities L/ic = 9, 
11, 13, 15, 20, 30, 40, 50, 60. The bending rigidity 
K is varied between 10~^ and 10~^ keeping a constant 
/i = 1. Successive boundary deformations are applied 
in logarithmically-spaced steps in order to explore the 
strain elastic response at both small and large strains. In 
Fig.Illis depicted a relatively dense network (L/ic = 15). 
Under a purely affine shear, the filaments initially ori- 
ented near 45° to the strain direction are stretched (in- 
dicated in red) and aligned, while those oriented near 
135° are compressed (blue) and eventually buckle. At 
the filament scale buckling occurs on wavelengths shorter 
than L due to the connectivity with the surrounding net- 
work through the cross-links, in contrast to classical Eu- 
ler buckling at the fundamental mode[T3]. 

In Fig. [2] are shown the normal and shear stresses, cta? 
and as, for a single network {L/£c = 13) in response 
to an imposed strain 7 for various bending stiffnesses k. 
The shear response starts linearly and eventually stiff- 
ens weakly, while the normal response is negative (the 
network wants to contract). For an ensemble-averaged 
network with isotropically oriented filaments, the nor- 
mal response must, by symmetry, be an even function 
of the shear deformation. Finite size effects, however, 
result in a small violation of this for any specific net- 
work studied. The shear stress increases with increasing 
K while the magnitude of the normal component exhibits 
the opposite trend. This behavior is quaUtatively consis- 
tent with the previous suggestion that the normal stresses 
arise from the asymmetry of the extensional response of 
single filaments [7] , since stiffer filaments are expected to 
bend less|T31 [TS] and the energy in Eq. (fTl) is symmet- 
ric. Increased bending, and eventual buckling, lead to an 
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FIG. 2: Raw shear (solid) and normal (dashed) stresses versus 
strain for different bending stiffnesses. 



increasingly asymmetric extensional response, with soft- 
ening under compression 1 1 6 j . 

In Fig. [3k we plot the shear modulus G = as/"/ and dif- 
ferential shear modulus K — das/dj normalized by the 



affine linear modulus Gaff — 






3 I , which 



corresponds to purely compression/extension deforma- 
tion without bending filaments [14|. This represents an 
upper bound for G and K at small strains, since bend- 
ing modes can only lower the energy for given bound- 
ary conditions. For these small strains, a linear regime 
with constant K = G = Gq is seen in all cases, while 
nonlinearities such as stiffening and/or softening occur 
at higher strains. The linear modulus Go increases both 
with density and bending stiffness in agreement with pre- 
vious work [m [15l [17] . Specifically, Go oc k for floppy 
filaments, while Gq saturates to Gaff for stiff filaments 
(inset) or high densities. 

In Fig. [4] we directly compare the shear and normal 
stresses, by plotting 0^ vs as, both of which have been 
normalized by Gaff- For stiffer filaments we see an ap- 
proximate quadratic dependence, which is to be expected 
by symmetry: since a;^ must be an even function of 
strain, it is generally expected to have an initial quadratic 
dependence for small strain, whereas as is expected to 
be approximately linear in strain. Interestingly, for more 
flexible filaments, this quadratic regime shrinks and an 
approximately linear dependence (— (Ttv — its '--^ 7 is ob- 
served, in which the two stresses are of comparable mag- 
nitude. This is consistent with what was found experi- 
mentally in Ref. [?]• This behavior can be understood in 
simple terms for highly asymmetric extensional response 
of filaments, in which both as and a^ are dominated by 
the stretched filaments oriented near 45° to the strain 
direction'?'. Our results suggest that the crossover from 
quadratic ctjv ~ a^ behavior occurs for smaller strains 
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FIG. 3: In (a), (b) and (c), the differential shear modulus 
K normalized by Gaff, the ratio of normal to shear stress, 
and the buckling order parameter B are all plotted versus 
shear stress for constant density L/£c — 15 and for various 



bending stiffnesses k = 

(•), K = 10-* (♦), « 
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10"^ (•). In the inset of (a), the 



ratio Go/Gaff is plotted versus the bending stiffness k for 
different densities L/lc = 9 (•), 13 (►), 20 (■) and 40 (♦). 
Here, Go grows linearly with k for soft filaments and saturates 
to Gaff for stiff ones; Go/Gaff also increases with density. 
Insets in (b) and (c) show data collapse for —onIos and B 
when plotted against os/i^- In (d) B is plotted for constant 
K, — 10^'' and various densities L/£c — 11 (*), L/£c = 13 (♦), 
L/4 = 15 (.), L/4 = 20 (►), L/4 = 30 (a), L/4 = 40 {<), 
L/lc = 50 (T) and L/£c = 60 (•). The inset in (d) shows the 
collapse of these data when plotted against as£c- Together 
with the insets in (b) and (c), this shows that the onset of 
nonlinearity coincides with the buckling transition. 



with more flexible filaments, which can be tested experi- 
mentally. 

To test this more directly, we plot the ratio —<Jn/<^s 
vs 7 in Fig. |5] Here, we see that the small-strain regime, 
characterized by um ^ 7^ and a^/crs ~ 7, decreases with 
increasing flexibility of filaments or decreasing network 
concentration. Interestingly, this is opposite to the pre- 
dicted behavior for thermally fluctuating networks in Ref. 
^. Thus, the normal stresses may provide an experimen- 
tal signature for thermal vs non-thermal systems. We 
also observe an apparent convergence of all networks for 
increasing strain, to a regime in which the two stress com- 
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FIG. 4: Normal stress Uiv versus shear stress as (both nor- 
malized by Gaff) for density L/£c = 13 and various filament 
bending stiffness. For stiff filaments, ctat -^ 7^ ^ a^ over 
a large initial stress region. For softer filaments this region 
becomes smaller so that on this linear scale it appears that 

CTJV ~ <JS- 
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FIG. 5: The ratio of normal to shear stress versus applied 
strain 7 for constant density L/£c — 15 and various filament 
bending stiffnesses. On decreasing k a peak grows, becom- 
ing more pronounced and moving to smaller strain. For large 
strain the curves depend weakly on k, showing a regime dom- 
inated by stretching only. The inset shows the same ratio 
— ctat/cts for constant stiffness k = 10-"" and different densi- 
ties. 



ponents are comparable. As the filaments become softer 
to bending, we observe a peak, with normal stresses ex- 
ceeding shear stresses for small strains. This may provide 
an explanation for the observation in Ref. ;2i, where the 
magnitude of normal stresses exceeded the shear stresses. 
As argued in Ref. [7], one possible origin of negative 
normal stresses can be the asymmetric extensional re- 



sponse, or force-extension relation of the constituent fil- 
aments. Even though the individual filaments in our 
model are assumed to be Hookean springs, with a sym- 
metric force extension curve, the fact that they can 
bend and buckle effectively gives rise to a softer re- 
sponse to compression and an asymmetric force-extension 
relationflGJ. To test whether this can explain our results, 
we measure the relative importance of buckling with an 
order parameter B defined by ^^cy'd^oioT jjgj-g^ tJig ^u- 
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merator refers to the incremental change in the elastic 
compression energy of the filament segments under com- 
pression in the network for a small strain step d^, evalu- 
ated as a function of the initial state of network strain 7. 
The denominator refers to the incremental change in to- 
tal elastic energy. This is normalized by its initial, small- 
strain value, where no buckling is expected. Thus, this 
order parameter should have initial value of unity, and 
should then decrease as filaments begin to buckle. We 
plot this in Fig.p] along with both K/Gaff and —(JnI<^s 
for several different networks of the same density, but 
with differing bending stiffness. We see that the onset of 
significant buckling, indeed, coincides with the onset of 
nonlinearities in the rheology. Furthermore, classical Eu- 
ler buckling theory predicts that the threshold force for 
buckling is proportional to k. By plotting both B and 
—ctn I^^S versus crgl n, we see collapse of the curves for 
various k (see insets), consistent with the dominant role 
of buckling in controlling the development of large nor- 
mal stresses. In Fig. [sli we plot B for various Ljtc and 
we see that the curves collapse (inset) by plotting against 
<Js^c^ confirming that the buckling order parameter actu- 
ally accounts for the collective buckling of intercrosslink 
segments [TB]. 

We have simulated networks of elastic rods with up to 
200,000 degrees of freedom, and for a range of strains 
up to 2. We find that coincident with the develop- 
ment of other nonlinear elastic properties, such networks 
very generally exhibit negative normal stresses of magni- 
tude up to and even exceeding the corresponding shear 
stresses. This highly unusual material property arises as 
a cooperative effect of filaments whose individual elastic 
behavior is linear. Qualitatively, the normal stress be- 
havior we observe can be understood, however, in terms 
of the nonlinear force extension that results from Eu- 
ler buckling of the filaments under compression [71 [TB] . 
We demonstrate this directly by measuring the degree 
of buckling in our networks. This can explain, in part, 
the apparent generality of negative normal stresses re- 
ported in recent experiments on a number of different stiff 
biopolymer systems [7]. Furthermore, our results suggest 
a possible way to distinguish experimentally between al- 
ternative theories of thermal [2l|3j [19] vs athermalfTZl ori- 
gins of nonlinear elasticity in biopolymer networks: we 
find that, the characteristic strain for the onset of nega- 
tive normal stress increases with network concentration, 
in contrast with the prediction for thermal networks [7] . 
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